Smooth-Particle Phase Stability with density and density-gradient potentials

Stable fluid and solid particle phases are essential to the simulation of continuum fluids and solids using Smooth Particle Applied Mechanics. We show that density-dependent potentials, such as Phi=(1/2)Sum (rho-rho_0)^2, along with their corresponding constitutive relations, provide a simple means for characterizing fluids and that a special stabilization potential, Phi=(1/2)Sum (delrho)^2, not only stabilizes crystalline solid phases (or meshes) but also provides a surface tension which is missing in the usual density-dependent-potential approach. We illustrate these ideas for two-dimensional square, triangular, and hexagonal lattices.Comment: 10 pages, 5 figure

Relativistic Hydrodynamic Flows Using Spatial and Temporal Adaptive Structured Mesh Refinement

Astrophysical relativistic flow problems require high resolution three-dimensional numerical simulations. In this paper, we describe a new parallel three-dimensional code for simulations of special relativistic hydrodynamics (SRHD) using both spatially and temporally structured adaptive mesh refinement (AMR). We used the method of lines to discretize the SRHD equations spatially and a total variation diminishing (TVD) Runge-Kutta scheme for time integration. For spatial reconstruction, we have implemented piecewise linear method (PLM), piecewise parabolic method (PPM), third order convex essentially non-oscillatory (CENO) and third and fifth order weighted essentially non-oscillatory (WENO) schemes. Flux is computed using either direct flux reconstruction or approximate Riemann solvers including HLL, modified Marquina flux, local Lax-Friedrichs flux formulas and HLLC. The AMR part of the code is built on top of the cosmological Eulerian AMR code {\sl enzo}. We discuss the coupling of the AMR framework with the relativistic solvers. Via various test problems, we emphasize the importance of resolution studies in relativistic flow simulations because extremely high resolution is required especially when shear flows are present in the problem. We also present the results of two 3d simulations of astrophysical jets: AGN jets and GRB jets. Resolution study of those two cases further highlights the need of high resolutions to calculate accurately relativistic flow problems.Comment: 14 pages, 23 figures. A section on 3D GRB jet simulation added. Accepted by ApJ

Tidal exposure or microhabitats: what determines sandy-beach nematode zonation? A case study of a macrotidal ridge-and-runnel sandy beach in Belgium

Lately, across-shore zonation has been found to be more important in structuring the nematode community of a tropical macrotidal sandy beach than microhabitat heterogeneity. To evaluate whether this zonation pattern applies to a temperate beach, a macrotidal ridge-and-runnels sandy beach in the North Sea was studied. We investigated whether a similar zonation occurs in sandbar and runnel microhabitats, and whether the runnels harbour a different community from the subtidal. Our results indicate that nematode communities from runnel and sandbar habitats are significantly different. In addition, horizontal zonation patterns for nematode communities differ between both habitats. Nematode assemblages from sandbars are divided to lower, middle and upper beach while upper and middle runnels cluster together. The subtidal and upper runnels showed dissimilar nematode assemblages, although runnels showed the same dominant species (Daptonema normandicum), which increases its abundance towards the upper runnels. This study illustrates the importance of microhabitat heterogeneity, which resulted in different zonation patterns across the sandy beach examined. The divergent zonation between sandbars and runnels in the macrotidal temperate sandy beach, compared with the pattern observed for a subtropical sandy beach with similar morphodynamics, indicates that generalizations about nematode distribution patterns should be made with caution

Solving One Dimensional Scalar Conservation Laws by Particle Management

We present a meshfree numerical solver for scalar conservation laws in one space dimension. Points representing the solution are moved according to their characteristic velocities. Particle interaction is resolved by purely local particle management. Since no global remeshing is required, shocks stay sharp and propagate at the correct speed, while rarefaction waves are created where appropriate. The method is TVD, entropy decreasing, exactly conservative, and has no numerical dissipation. Difficulties involving transonic points do not occur, however inflection points of the flux function pose a slight challenge, which can be overcome by a special treatment. Away from shocks the method is second order accurate, while shocks are resolved with first order accuracy. A postprocessing step can recover the second order accuracy. The method is compared to CLAWPACK in test cases and is found to yield an increase in accuracy for comparable resolutions.Comment: 15 pages, 6 figures. Submitted to proceedings of the Fourth International Workshop Meshfree Methods for Partial Differential Equation

The Chemical Compositions of the Type II Cepheids -- The BL Her and W Vir Variables

Abundance analyses from high-resolution optical spectra are presented for 19 Type II Cepheids in the Galactic field. The sample includes both short-period (BL Her) and long-period (W Vir) stars. This is the first extensive abundance analysis of these variables. The C, N, and O abundances with similar spreads for the BL Her and W Vir show evidence for an atmosphere contaminated with $3\alpha$-process and CN-cycling products. A notable anomaly of the BL Her stars is an overabundance of Na by a factor of about five relative to their presumed initial abundances. This overabundance is not seen in the W Vir stars. The abundance anomalies running from mild to extreme in W Vir stars but not seen in the BL Her stars are attributed to dust-gas separation that provides an atmosphere deficient in elements of high condensation temperature, notably Al, Ca, Sc, Ti, and $s$-process elements. Such anomalies have previously been seen among RV Tau stars which represent a long-period extension of the variability enjoyed by the Type II Cepheids. Comments are offered on how the contrasting abundance anomalies of BL Her and W Vir stars may be explained in terms of the stars' evolution from the blue horizontal branch.Comment: 41 pages including 11 figures and 4 tables; Accepted for publication in Ap

Effects of a Supermassive Black Hole Binary on a Nuclear Gas Disk

We study influence of a galactic central supermassive black hole (SMBH) binary on gas dynamics and star formation activity in a nuclear gas disk by making three-dimensional Tree+SPH simulations. Due to orbital motions of SMBHs, there are various resonances between gas motion and the SMBH binary motion. We have shown that these resonances create some characteristic structures of gas in the nuclear gas disk, for examples, gas elongated or filament structures, formation of gaseous spiral arms, and small gas disks around SMBHs. In these gaseous dense regions, active star formations are induced. As the result, many star burst regions are formed in the nuclear region.Comment: 19 pages, 11 figures, accepted for publication in Ap

SPH Simulations of Counterrotating Disk Formation in Spiral Galaxies

We present the results of Smoothed Particle Hydrodynamics (SPH) simulations of the formation of a massive counterrotating disk in a spiral galaxy. The current study revisits and extends (with SPH) previous work carried out with sticky particle gas dynamics, in which adiabatic gas infall and a retrograde gas-rich dwarf merger were tested as the two most likely processes for producing such a counterrotating disk. We report on experiments with a cold primary similar to our Galaxy, as well as a hot, compact primary modeled after NGC 4138. We have also conducted numerical experiments with varying amounts of prograde gas in the primary disk, and an alternative infall model (a spherical shell with retrograde angular momentum). The structure of the resulting counterrotating disks is dramatically different with SPH. The disks we produce are considerably thinner than the primary disks and those produced with sticky particles. The time-scales for counterrotating disk formation are shorter with SPH because the gas loses kinetic energy and angular momentum more rapidly. Spiral structure is evident in most of the disks, but an exponential radial profile is not a natural byproduct of these processes. The infalling gas shells that we tested produce counterrotating bulges and rings rather than disks. The presence of a considerable amount of preexisting prograde gas in the primary causes, at least in the absence of star formation, a rapid inflow of gas to the center and a subsequent hole in the counterrotating disk. In general, our SPH experiments yield stronger evidence to suggest that the accretion of massive counterrotating disks drives the evolution of the host galaxies towards earlier (S0/Sa) Hubble types.Comment: To appear in ApJ. 20 pages LaTex 2-column with 3 tables, 23 figures (GIF) available at this site. Complete gzipped postscript preprint with embedded figures available from http://tarkus.pha.jhu.edu/~thakar/cr3.html (3 Mb

The mass function

We present the mass functions for different mass estimators for a range of cosmological models. We pay particular attention to how universal the mass function is, and how it depends on the cosmology, halo identification and mass estimator chosen. We investigate quantitatively how well we can relate observed masses to theoretical mass functions.Comment: 14 pages, 12 figures, to appear in ApJ

The Structure of Isothermal, Self-gravitating Gas Spheres for Softened Gravity

A theory for the structure of isothermal, self-gravitating gas spheres in pressure equilibrium in a softened gravitational field is developed. The one parameter spline softening proposed by Hernquist & Katz (1989) is used. We show that the addition of this extra scale parameter implies that the set of equilibrium solutions constitute a one-parameter family, rather than the one and only one isothermal sphere solution for Newtonian gravity. We demonstrate the perhaps somewhat surprising result that for any finite choice of softening length and temperature, it is possible to deposit an arbitrarily large mass of gas in pressure equilibrium and with a non-singular density distribution inside of r_0 for any r_0 > 0. The theoretical predictions of our models are compared with the properties of the small, massive, quasi-isothermal gas clumps which typically form in numerical Tree-SPH simulations of 'passive' galaxy formation of Milky Way sized galaxies. We find reasonable agreement despite the neglect of rotational support in the models. We comment on whether the hydrodynamical resolution in our numerical simulation of galaxy formation is sufficient, and finally we conclude that one should be cautious, when comparing results of numerical simulations involving gravitational softening and hydrodynamical smoothing, with reality.Comment: 22 pages Latex + 12 figure

Tidal disruption of dark matter halos around proto-globular clusters

Tidal disruption of dark matter halos around proto-globular clusters in a halo of a small galaxy is studied in the context of the hierarchical clustering scenario by using semi-cosmological N-body/SPH simulations assuming the standard cold dark matter model ($\Omega_0 = 1$). Our analysis on formation and evolution of the galaxy and its substructures archives until $z = 2.0$. In such a high-redshift universe, the Einstein-de Sitter universe is still a good approximation for a recently favored $\Lambda$-dominated universe, and then our results does not depend on the choice of cosmology. In order to resolve small gravitationally-bound clumps around galaxies and consider radiative cooling below $T = 10^4 K$, we adopt a fine mass resolution (m_{\rm SPH} = 1.12 \times 10^3 \Msun). Because of the cooling, each clump immediately forms a `core-halo' structure which consists of a baryonic core and a dark matter halo. The tidal force from the host galaxy mainly strips the dark matter halo from clumps and, as a result, theses clumps get dominated by baryons. Once a clump is captured by the host halo, its mass drastically decreases each pericenter passage. At $z = 2$, more than half of the clumps become baryon dominated systems (baryon mass/total mass $> 0.5$). Our results support the tidal evolution scenario of the formation of globular clusters and baryon dominated dwarf galaxies in the context of the cold dark matter universe.Comment: 9page, 13 figures. Accepted for publication in ApJ. A high-resolution PDF of the paper can be obtained from http://th.nao.ac.jp/~takayuki/ApJ05